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Low-temperature neutron diffraction and NMR studies of field-induced phases in linarite are 
presented for magnetic fields H || b axis. A two-step spin-flop transition is observed as well as 
a transition transforming a helical magnetic ground state into an unusual magnetic phase with 
sine-wave modulated moments || H. An effective J 1 -J 2 single-chain model with a magnetization- 
dependent frustration ratio a e s = —J 2 /J 1 is proposed. The latter is governed by skew interchain 
couplings and shifted to the vicinity of the ferromagnetic critical point. It explains qualitatively the 
observation of a rich variety of exotic (for strongly correlated cuprate spin-1/2 Heisenberg systems) 
longitudinal collinear spin-density wave SDW P states (9 > p > 2). 

PACS numbers: 75.10.Jm, 75.25.-j, 75.30Kz, 75.40Cx 


Recently, frustrated spin chains with ferromagnetic 
nearest neighbor J\ (FM-NN) and antiferromagnetic 
2 nd neighbor J 2 (AFM-NNN) exchange have been dis¬ 
cussed in the context of novel states of matter. Close 
to the saturation field by tuning the frustration ratio 
a = —J 2 /J 1 of an isotropic model, a sequence of dis¬ 
tinct spin-multipolar (MP) phases should develop with 
well-defined phase boundaries, which are described as a 
Tomonaga-Luttinger-liquid of p-magnon bound states [1 
9]. They compete with exotic longitudinal spin density 
wave (SDWp) correlations, which should prevail in lower 
magnetic fields. Interchain exchange should weaken the 
MP correlations, while magnetic anisotropy should sta¬ 
bilize them [6-9]. 

A proof of existence for spin-MP ordering in real quasi- 
ID materials is still lacking. The FM-AFM chain J 1 -J 2 
compound LiCuV 04 was considered as a candidate, un¬ 
dergoing a transition into an incommensurate (ICM) he¬ 
lical phase below 2.1 K, into a spin-flop phase in fields of 
~2.5T and an exotic SDW 2 phase above ~8T [10-13]. 
A shift of the SDW 2 propagation vector, consistent with 
longitudinal density waves of bound p = 2-magnons, was 
reported for low fields of 8-14.5 T, together with a transi¬ 
tion from long- to short-range magnetic order [13]. It was 
interpreted as a signature of coexisting SDW 2 and bond- 
nematic order, a view disputed in Refs. [9, 14] and instead 
related to a pinned SDW 2 . Via magnetization and NMR 
it was concluded that MP correlations in LiCuVCH can 
exist only in a narrow high field range ~ 40 T [5, 15, 16]. 

The issue not resolved in this dispute is the relation¬ 
ship between SDW and MP in the long-range ordered 
phases appearing in 2D and 3D (in ID the precursor 
’’phases” overlap [3]). Theoretically, the possibility of ho¬ 


mogeneously coexisting SDW 2 and nematic phase and/or 
phase separation has been suggested for the isotropic 
model in an extreme quasi-ID regime for specific intra¬ 
chain and very weak 2D FM interchain couplings based 
on perturbative scattering theory [17]. In contrast, only 
a l st -order phase transition was predicted for the same 
model [9,14] . Also predicted is a phase separation and/or 
a l st -order phase transition between nematic and FM 
phase in a 3D bcc structure for the same J 1 -J 2 model 
and a similar approach [18] as in Ref. [17]. The situation 
is far from clear in LiCuVCb. The phase diagram has 
not been studied in detail due to the high fields required 
to access it [19]. Further, the influence of defects on the 
magnetic properties is not well understood [16, 20, 21]. 
Thus, what is lacking in this context is a comprehensive 
study of a clean frustrated FM-AFM spin chain material 
to properly define these issues. 

A unique example of a frustrated spin chain system for 
such studies is linarite. It crystallizes in the monoclinic 
space group P2i/m [22], forming buckled CuC >2 chains 
along the b axis. These have been modelled as a s = | 
spin-chain with FM-NN Ji = —100 K and AFM-NNN 
J 2 = 36 K [23]. In this J-parameter range the saturation 
field is about 10 T, allowing full experimental access to 
the magnetic phase diagram. For a magnetic field H || b 
axis the magnetic phase diagram contains five different 
regions I (elliptical helix) to V [24-26] (see Supplement). 
Region V displays very weak thermodynamic signatures, 
and it was unclear, whether it is a distinct thermody¬ 
namic phase. 

Here, we fully characterize its field-induced phases by 
means of neutron diffraction (ND) and 4 H nuclear mag¬ 
netic resonance (NMR). We establish the magnetic or- 
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dering vectors, and that region V represents a thermody¬ 
namic phase. For phase V, we determine the field depen¬ 
dence of the ICM SDW ordering vector and newly dis¬ 
cover complex states which might be understood in terms 
of phase separation between MP and SDW P states. 

ND was carried out using the single crystal instrument 
DIO at the Institute Laue Langevin, France, and the in¬ 
strument Wombat at ANSTO, Australia. For the DIO ex¬ 
periment, the sample from a previous study was used [24]. 
A second single crystal of linarite (9 x 3 x 0.5 mm 3 ) from 
the Grand Reef Mine, Arizona, was used for the experi¬ 
ment on Wombat. The samples were placed in cryomag- 
nets with maximum field/base temperature of 6T/1.7K 
(DIO) and 12T/1.5K (Wombat), with the magnetic field 
applied along the crystallographic b axis. With this setup 
and a neutron wavelength of 2.36 A we were restricted 
along the b direction to —0.25 < k < 0.25 in reciprocal 
space for the DIO experiment, while for the experiment 
on Wombat a wavelength of 4.61 A was used resulting in 
an access range —0.19 < k < 0.19. 

1 H-NMR ( 3 y = 42.5749 MHz/T) studies were per¬ 
formed for T < 2.8 K using a phase-coherent Tecmag 
spectrometer in combination with a He-flow cryostat. 
Frequency scans were conducted down to 1.7 K and at 
external fields H || b between 1.5 T and 7.5 T. The same 
single crystal was used as for the D10 ND study. All 
NMR spectra were collected using a n/2 — t — tt Hahn 
spin-echo pulse sequence. The spectra have not been cor¬ 
rected by the tiny spin-spin relaxation time T 2 ~ 10/zs. 

By ND, scans of (0 k 0.5) with varying k were per¬ 
formed at 1.7 K for fields up to 6T. For this temperature 
the sequence of phases I III -IV is traversed with increas¬ 
ing held. Consistent with the 0T propagation vector, 
the magnetic Bragg peak (0 —0.186 0.5) is observed at 
low fields (Fig. la). For increasing magnetic held, at 
the boundary to phase III a second commensurate (CM) 
magnetic Bragg peak appears at (0 0 0.5), corresponding 
to spins coupled parallel along the a and b axes and an¬ 
tiparallel along c. In the held range 2.60 to 2.95 T both 
peaks coexist, while for higher helds (in phase IV) only 
the commensurate Bragg peak remains. 

The coexistence of two spin structures in phase III is 
also observed in NMR. Fig. lb displays the absolute res¬ 
onance change v ~ v 0 of the spectra taken in phase I, 
III, and IV (for details see Supplement). The six NMR 
lines in phase I are understood in terms of an ICM helical 
structure [26]. In phase IV four discrete resonance lines 
were observed. For linarite the only CM spin structure 
with four resonance lines is an AFM, propagating along 
the a or c direction, consistent with ND. Conversely, in 
phase III in total 12 resonance lines are detected; four 
with high intensity at small shifts v — fq = —1.4 to 
0.5 MHz, eight lines for shifts larger than ±1.5 MHz. The 
four central lines resemble those of phase IV. The remain¬ 
ing eight peaks result from a modification of the helical 
phase I structure. Thus, a phase separation occurs into 


two spin structures. The two phases compete, as for in¬ 
creasing field the intensity of the central lines increase 
while the other decrease, reflecting a growth of the phase 
volume of the first on behalf of the second. The same 
behavior is seen for the field dependence of the corre¬ 
sponding Bragg peak intensities in ND (Fig. la). 

To determine the magnetic structure with propagation 
vector k = (0 0 0.5) of phase IV, the intensity of 33 mag¬ 
netic Bragg peaks (20 inequivalent) was measured at 4T 
by ND. A refinement of the data (Rf = 9.5 %) using the 
program FullProf [27] reveals that the spins are lying in 
the ac plane, with an angle of —27° off the a axis (roughly 
parallel to [101]), the same as one of the spin components 
of the phase I helix [24] . From the refinement an ordered 
moment of 0.79(1)/xb per Cu atom is derived. Similar re¬ 
finements at 5.5 T (20 inequivalent Bragg peaks) yield the 
same spin structure with a moment of 0.73(2)/zb per Cu 
atom (Rp = 12.5%). The decrease of the AFM moment 
with held, and the observation of small field-induced FM 
contributions on top of nuclear Bragg peaks, reflects the 
development of field-induced spin polarization. 

For the determination of spin structures in phase III, 
two sets of magnetic Bragg peaks ( hkl)M were collected 
at 2.8 T using the relation ( hkV)M — (hkl)pr ± k. For the 
ICM structure the propagation vector k = (0 0.186 0.5) 
was used, for the CM structure k = (0 0 0.5). The CM 
structure in phase III is refined using 15 peaks (14 in¬ 
equivalent) with the phase IV spin model. The refine¬ 
ment of 18 inequivalent Bragg peaks of the ICM structure 
yields a circular helix structure (Rp = 14.5 %), where the 
moments of 0.64(2)/xb lie roughly in the be plane. 

In the related chain systems LiCuV 04 and LiCu 2 02, 
applying magnetic held rotates the normal of the helical 
structure parallel to the held. Here, such a spin hop of 



FIG. 1. a.) Neutron scattering scans for linarite along k at 
1.7 K at different fields, reflecting the crossing from phase I 
via III into IV. Solid lines are Gaussians fitted to the data 
to determine peak positions, b.) NMR spectra for linarite at 
1.7K at different fields (spectra offset for clarity). 
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FIG. 2. a.) The magnetic Bragg peak position (0 k 0.5) shift¬ 
ing with field in phase V. b.) Field dependence of the prop¬ 
agation vector in phase V compared to the ID model SDW P 
states with 2 < p < 9. Theory lines include the magnetization 
obtained from M(H,T ) scans [25]. 


the helix into the ac plane is prohibited by the mono¬ 
clinic angle /3. Instead the spins flop into a collinear spin 
structure in the ac plane. In phase III the spins start to 
flop into the ac plane forming a collinear spin arrange¬ 
ment, while a coexisting helical phase is retained. The 
fact that in phase III a circular helix (spinning plane in 
the be plane) replaces the elliptical helix reflects that for 
the latter it is energetically costly to keep the large mo¬ 
ment axis aligned along the field direction. For larger 
fields all spins are flopped into the ac plane forming the 
collinear phase IV. 

Next, we have performed ND via k scans through re¬ 
gion V (Fig. 2a). Surprisingly, and in spite of the very 
weak signatures defining this region in thermodynamic 
measurements [24, 25] , we observe magnetic Bragg peaks 
of the same width as the nuclear peaks. These peaks 
are even observed in the intermediate field regime ^4T, 
where in thermodynamic measurements no anomalies 
were detected [25]. It implies that region V is a distinct 
and highly unusual thermodynamic phase. For magnetic 
structure determination, a set of 8 inequivalent Bragg 
peaks was collected at 6T. A refinement yields a sine- 
wave modulated structure, with the spins aligned parallel 
to the b axis ( Rf = 7%). Surprisingly, the SDW mag¬ 
netic moment amplitude is only 0.44(1)^b- This value is 
much smaller than what would be expected from an ex¬ 
trapolation of the field dependence of the magnetic mo¬ 
ments measured in phase I and IV. Further, the propa¬ 
gation vector (0 —k 0.5) shifts in k with field, as shown 
in Fig. 2. 


The sine-wave modulated structure with moments 
along the field direction together with the shift of the 
k value reminds of the prediction of the longitudinal 
collinear SDW within hard-core boson approximation 
[1, 3, 8], where the shift depends on the number of bound 
magnons p in the coexisting or neighboring MP phase: 

kyd (1 — M/Ms) ' 

7T p 

Here, d denotes the distance of neighboring Cu spins 
along the b axis, Ms is the saturation magnetization 
[28, 29]. To compare Eq. 1 with the situation for linarite, 
the curves with various p values are included in Fig. 2b. 
Surprisingly, at first glance no agreement is found over 
a wide field range between the experimentally observed 
evolution of k y and the theoretical prediction for a single 
chain with fixed field independent exchange interactions 
for any fixed value of p. However, taking into account the 
interchain coupling we will arrive at a magnetization de¬ 
pendent reduced effective value a e s{M/M s ). Then, the 
different p values seen experimentally can be qualitatively 
understood (see below and Supplement). 

Static magnetic order in phase V of linarite is also 
observed by NMR. 1 H-NMR frequency scans were per¬ 
formed in the field range 3-7.5T with an increment of 
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FIG. 3. The 1 H-NMR spectra of linarite for different magnetic 
fields and temperatures in phase V and in the paramagnetic 
polarized state at 7.5 T close to saturation. 
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0.5 T at different temperatures. At 2.8 K a paramagnetic 
signal is observed which is composed of two almost over¬ 
lapping lines from two inequivalent *H sites [26]. For 
fields < 6 T, upon lowering T below a critical value Ty 
the spectrum develops horn-shaped NMR characteristics, 
that is two distinct peaks with a finite intensity in be¬ 
tween (Fig. 3a-c). It can be accounted for by the SDW 
structure with only a magnetic component along the b 
direction (compare LiCuV 04 [16]). The transition tem¬ 
peratures derived from NMR match those of phase V 
obtained from thermodynamics for H || b axis [24, 25], 
and define the phase boundary in the field range 3.5 -6 T, 
where no transition has been detected in thermodynamic 
quantities. Our findings imply that phase V encloses all 
other magnetic phases (see Supplement). 

Increasing the field to above 6T within phase V (in 
Fig. 3d/e: 1.7 K) produces a transfer of spectral weight 
from the horn-shaped structure to a broadened two-peak 
structure appearing in the middle of the SDW pattern. 
The shift of the latter with increasing field follows the 
shift of the paramagnetic polarized NMR signal close to 
saturation. Qualitatively, this implies the presence of two 
different local environments in the sample, viz., a phase 
separation occurs. In part of the sample there is SDW or¬ 
dering producing the horn-shaped NMR spectra. In con¬ 
trast, the regions of the sample exhibiting the broadened 
two-peak structure show no static magnetic order. Ac¬ 
cording to Refs. [1-3] for coupled frustrated spin chains 
a field-induced transition from the SDW P into a p-MP 
phase is expected in high fields. This transition should 
appear as one from a magnetically long-range ordered 
into one without static dipolar long-range order. Hence, 
we suggest that the phase separation observed in phase 
V is related to the transition from a SDW p phase into 
one with dominant MP character in a quasi-ID material. 

To discuss our experimental results, we employ weakly- 
coupled J\ — Jn. chains in a magnetic field H along the z 
axis. The Hamiltonian reads as 

H = J\ ^ S M • Sz,i+i + J 2 S Zji • S/ ii+2 

l,i l,i 

+ HY,S* + J ic J2 S M-S( 2 ) 

where S i t i is a spin-i operator at site i in chain l and J lc is 
a diagonal interchain coupling (see Fig. 4(a)). As shown 
above, the ICM propagation vector along the chain is 
k y = 0.1867T at H = 0; however, a single Ji — J 2 chain 
with a = J 2 /IJ 1 I = 0.36 gives k y « 0.36 77T. This dis¬ 
crepancy can be resolved by taking a specific diagonal 
J ic « 10K = 0.11 Ji|. Theoretically, the propagation vec¬ 
tor in the single J 1 -J 2 chain is found from the maximum 
position of the static spin-structure factor S(k). Due to 
strong quantum fluctuations and the resulting absence of 
static magnetic order in ID with short range-couplings, 


only, it cannot be found from < Sf >. But as a pre¬ 
cursor, the former reflects nevertheless the SDW mod¬ 
ulations (induced by the coupling to neighboring chains 
in 2D and 3D) we are looking for here. This maximum 
position of S(k) is reduced by decreasing a and it ap¬ 
proaches 0 at the FM critical point a c = 1/4. Such a 
reduction of the propagation vector is also realized by 
increasing J; c at fixed a for a system of coupled chains 
(see Fig. 4a and Fig. 5 in Ref. [30]). It is thus inter¬ 
preted that the interchain coupling reduces the value of 
a. Since a single chain is more convenient than cou¬ 
pled chains for the discussion of MP states, we first per¬ 
formed a mapping from two coupled J 1 -J 2 chains with 
a = 0.36 and J lc = 0.1|Ji| with periodic perpendicular 
boundary conditions onto an effective single J 1 -J 2 chain 
with a e ff = —J 2 /J 1 = a e s(a, M/M s ). For a wide range 
of the magnetization M/M s the values of a e ff were es¬ 
timated by fitting the dynamical spin-structure factors 
S(q, uj). An example is shown in Fig. 4(b) (for more de¬ 
tails see Supplement). The estimated values of a e (t are 
plotted vs. M/M s in Fig. 4(c). Note the vicinity to a c . 

Next, we found the number of bound magnons p for a 
given a by calculating the binding energy of a p-magnon 
bound state near the saturation field, which is defined as 



site index/ 


site index/ 


FIG. 4. (a) Lattice model of weakly-coupled Ji — J 2 chains, 
(b) Fitting of S(q,ui) at M/M s = 0.2; between q = 
for coupled two-chain with a = 0.36, Jic = 0.1|Ji| and q = n 
for a single chain with a e ff = 0.257. (c) Phase diagram of 
the multi-magnon bound states and the estimated effective 
frustration ratio a e fi for linarite. (d) Suggested propagation 
vector using Fig. 4(c). Local spin densities (Sf) at a e g = 0.26 
and M/Ms = 0.95 for (e) periodic and (f) open chains. 
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Eh(p) = ^ [E(S Z = S m , ax - p) - E(S Z = S max )] 

- [E(S Z = 5 max - 1) - £(S* = S max )], (3) 

where E(S Z = S) is the ground-state energy with the 2 ; 
component of the total spin S z = S , and S x = S max cor¬ 
responds to the fully polarized state. When the largest 
value of Eh(p){> 0) is given by p = p max , we can prove 
that the Pmax-magnon bound state is the most stable 
state; whereas, if £b(p) < 0 for all p values no low-energy 
magnon bound state exists. The results are shown in 
Fig. 4(c). According to Refs. [1-3], the value of p in¬ 
creases with approaching a e ff = 1/4 and the region of 
a e s becomes narrower for larger p. Based on the rela¬ 
tion between p and a e ff, the propagation vector of linar- 
ite is suggested to evolve similarly as shown in Fig. 4(d) 
based on calculations within our effective ID model, in 
semi-quantitative agreement with the experimental data 
(Fig. 2(b)). For a brief discussion of the last right dashed 
part of the broken red lines related to p = 2 obtained 
within an analogous xyz-anisotropic Heisenberg model 
(to be published elsewhere) see Supplement. 

The phase separation observed at high field may be at¬ 
tributed to the low effective frustration ratio a e ff ~ 1/4, 
i.e., close to the critical point. This means that the FM 
state is almost degenerate to other lower spin states. As 
an illustration, the local spin densities ( S?) with periodic 
and open boundaries at M/M s = 0.95 for a e ff = 0.26 are 
plotted in Figs. 4(e) and (f), respectively. A uniform 
distribution is naturally expected for periodic chains; 
whereas, interestingly, a phase separation into partial po¬ 
larized and unpolarized phases occurs by open chain ends 
or possibly also near strong enough impurities in the bulk 
being a common small disorder in real materials. 

To conclude, linarite exhibits a field-induced behav¬ 
ior generic for a FM-NN/AFM-NNN frustrated chain 
system. In addition, a two-step spin-flop transition is 
present for external magnetic fields applied along the 
b axis. Further, a longitudinal sine-wave modulated 
spin-structure phase encloses the other ordered phases. 
Here, at relative low fields above the helical phase a 
shift in the propagation vector qualitatively similar to 
LiVCuCH (with SDW 2 states, only) is observed. How¬ 
ever, to the best of our knowledge we report the first ob¬ 
servation of several exotic (for Heisenberg spin-1/2 sys¬ 
tems) collinear longitudinal SDW P states (by changing 
the external field), reaching even p = 9. We believe 
that this result is related to the appropriate bare value 
of a and the strong enough skew interchain coupling. 
Altogether, linarite appears to be a good candidate to 
show MP behavior. A more detailed and comprehensive 
study, its exotic SDW P states, and their interplay with 
field-induced phase separation and exchange anisotropy 
provides a challenge for future work. Also the compari¬ 
son with other rare cases of field-induced phase separa¬ 


tion in frustrated zigzag chain magnetic spiral systems 
as MnWCH [31] is of interest for deeper insights in their 
complex physics including multiferroicity. 
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In this Supplementary part we present the magnetic phase diagram of linarite PbCuS 04 ( 0 H )2 with H || b axis 
for the convenience of the reader once again, but now extended by the inclusion of new NMR data and under 
consideration of additional neutron scattering data taken in magnetic fields. These new data establish phase V as a 
distinct thermodynamic phase. Furthermore, the evolution of the 4 H-NMR spectra in phase III is shown as a function 
of the magnetic field. The nature of the different subspectra are discussed and compared to the 4 H-NMR spectra in 
the phases I and IV. Finally, we explain in more detail the mapping procedure of two coupled J\ — Ji, Ji C -chains onto 
an effective single J 1 -J 2 chain but with a magnetization dependent frustration ratio a e ff = Q? e ff (ct [ M/M s ]). 


THE MAGNETIC PHASE DIAGRAM OF LINARITE FOR H || b 

The phase diagram for a magnetic field H applied parallel to the crystallographic h axis of linarite contains five 
different regions/phases I to V [1-3] (see Fig. SI). The ground state phase I below ~ 2.5T is formed by an incom¬ 
mensurate (ICM) elliptical helix (ordering wave vector k = (0 0.186 0.5), moment size /j, or d ~ 0.64 to 0.83/is). The 
field-induced transition from phase I into phase IV involves either a hysteretic transition (region II, below ~ 0.6K), 
a direct one (up to ~ 1.25K), or the transit through a phase III (up to ~ 1.9K and 3.2T). Region V (up 9.5T) 
displays only very weak thermodynamic signatures, therefore it was not clear whether it is a distinct thermodynamic 
phase. This important issue has been settled by 4 H-NMR and neutron diffraction measurements presented in this 
manuscript. The transition temperatures into phase V obtained from the temperature dependence of our 1 H-NMR 
spectra are included in Fig. SI for the first time. 
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FIG. SI. Phase diagram of linarite with the applied field along the crystallographic b axis. For a brief description of the phases 
I to V see text and Refs. 1-3. Transition temperatures obtained from new 1 H-NMR data are depicted as red diamonds. 
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NATURE OF THE 1 H-NMR SPECTRA IN PHASE III 

Phase III is wedged in between phase I and phase IV in the magnetic phase diagram of linarite for H || b. From 
the thermodynamic studies it was speculated that phase III is part of a two-step spin spiral reorientation process by 
going from phase I via phase III to phase IV. This notion needs to be verified on a microscopic level by means of a 
combined NMR and neutron scattering study. Here, additional information on the NMR part is presented. 

1 H-NMR was measured at 1.7K in 2.5 and 3.0T and is shown in Fig. S2(a) and (b). For both fields 12 resonance 
lines are detected, four of them with a high intensity in the middle of the spectra close to the absolute resonance 
change v — z'l = 0 MHz, and eight lines situated at the borders of the spectra. The four high-intensity lines in the 
middle of the spectra closely resemble the shape of the spectra in phase IV, with the latter shown in Fig. S2(c) for an 
external field of 3.5 T. In turn, by increasing the external magnetic field from 2.5 to 3T, the whole spectrum shifts to 
lower frequencies v — ul- 

Remarkably, the four lines in the middle gain intensity with increasing field, whereas the remaining eight peaks lose 
intensity. At 2.5 T the four peaks have a relative intensity of 76%, which is increased to 96% at 3.0 T, as determined 
by integrating the subspectra and assuming constant T% times in 2.5 and 3T for the individual peaks. 

The eight-peak subspectrum appears to be discrete, which would point towards an underlying commensurate mag¬ 
netic structure. However, also here the spin-spin correlation time X 2 is very short, which would hamper identification 
of a horn-shape spectrum expected for an incommensurate spiral structure (cf. the situation for phase I [3]). Our neu¬ 
tron scattering data in phase III provide evidence for a coexistence of the phase IV commensurate magnetic structure 
with an incommensurate spin spiral. Then, on the supposition of horn-shaped spectra with pairs of lines reaching 
from one to the other border of the spectrum, viz., from negative to positive v — ^l, the intensity in between the peaks 
is easily suppressed due to T 2 being of the order of 10 n s, resulting in the appearance of seemingly discrete peaks. 

Under the assumption that the eight remaining peaks are related to the incommensurate helical ground state (due 
to the same number of resonance lines), a modification of this structure may be expected at the phase boundary I-III. 
The abrupt change observed in v — at the phase boundary from phase I to III makes a change of the spin structure 
mandatory. Modifications could be a significant change in (either) the spinning plane, the magnetic moment, or the 
propagation vector. Hence, the NMR data is fully consistent with the conclusion from neutron scattering, i.e., a 
coexistence in the sample of a spin spiral with phase IV regions. 
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FIG. S2. (a) and (b) 1 H-NMR spectra of linarite in phase III at T = 1.7 K for poH = 2.5 and 3.0 T || b axis: Four discrete lines 
with high intensities (arbitrarily labeled by red numbers) are visible in the middle of each spectrum, which can be associated 
to the spectra known from phase IV. Eight additional resonance lines (arbitrarily numbered in black) are placed at the border 
of each spectrum, belonging to a new magnetic structure but which is closely related to the helical magnetic ground state, (c) 
For comparison, 1 H-NMR spectrum of linarite in phase IV at T = 1.7 K for poH = 3.5 T || b axis. 
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MAPPING OF TWO COUPLED Ji - J 2 , J ic CHAINS ONTO AN EFFECTIVE SINGLE Ji - J 2 CHAIN 

In the main text, we have performed a mapping from the J\ — J 2 chains weakly coupled by the diagonal interchain 
exchange interaction J lc (see Fig. S3(a)), 

H = Jr ^ S M ■ S M+1 + J2^S lti -S l>i+2 + H ^ s f ii + J ic ^ S M -S Py , (SI) 

l,i l,i l,i l,V ,i,i' 

onto an effective single J\ — J 2 chain (see Fig. S3(b)), 

H = h Y, Si • Sj +1 +JiJ2 s >- S,;+2 + H Y, Si , (S2) 

i i i 

where S/,» is a spin -tj operator at site i in chain l. The objective of this mapping is to mimic the internal magnetic 
properties of a chain in the coupled J\ — J 2 model by an effective single J\ — J 2 chain. A reliable mapping can be 
achieved by calculating the dynamical spin structure factor 

S(q,w) = \{^\S(cO\ihMu - (K - E 0 )), (S3) 

v>0 

where the static structure factor 5(q) is simply the Fourier transform of Sjb, E v and |i 1> V ) are the z/-th eigenenergy 
and eigenstate of the system, respectively (the ground state denoted by v = 0). For some fixed magnetizations M/M s 
we determined the value of a e g = J 2 /\J\\ to reproduce the spin structure factor of the 15 x 2 coupled chains by 
the effective 15-site chain as much as possible. As an example of the mapping, the comparison of the static and 
dynamical structure factors between the original coupled chains and the effective single chain at M/M s = 0.6 is shown 
in Fig. S3(c) and (d). We can see a good agreement. 



FIG. S3. Lattice models of (a) weakly-coupled Ji — J 2 chains and (b) effective Ji — J 2 chain, (c), (d) Comparison of the static 
and dynamical spin structure factors between the weakly-coupled Ji — J 2 chains with a = J 2 /|Ji| = 0.36, Ji c /|A| = 0.1 (blue) 
and the effective single Ji — J 2 chain with a e s = J 2 /|Ji| = 0.263 (red) for a fixed magnetization M/M s = 0.6. The static 
structure factor of the coupled chains S(q x ) has been averaged over q y — 0 and q y = n. The orange and green dotted lines 
denote the lower and upper bounds of the spinon continuum, respectively. 
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Notice that the obtained ’’renormalized” a e g depends naturally on both the bare value a and on the type and 
strength of the interchain coupling. The smaller the former and the larger the skew interchain couplings are, the 
smaller is the resulting a e ff. To illustrate this point, we mention that for larger bare a values above 0.7, as relevant 
for LiVCu 04 , one remains in the p = 2-region for all magnetic fields in accord with the data of Mourigal et al. [4]. 
In contrast, for Li 2 Cu 02 with a smaller bare a ~ 0.32 as compared to linarite one reaches already the ferromagnetic 
phase with ct e ff < 0.25 [5] in accord with the observed alignment of magnetic momenta along the chains and without 
any anomalous longitudinal collinear SDW P states. Finally, since the pitch is only weakly affected by the exchange 
anisotropy, the inclusion of the last ’’zigzag” into the curve shown in Fig. 4(d) of the main text obtained from an 
anisotropic generalization of our single chain Hamiltonian is reasonable (see also the next subsection of the present 
Supplement). 


NEMATIC (p = 2) STATE NEAR THE SATURATION 


According to preliminary (unpublished) measurements for our title compound the material is highly anisotropic (see 
also Ref. [ 6 ]) and a more general spin model might be provided by a Heisenberg chain with xyz-exchange anisotropy 
for any a e ff- Then the Hamiltonian is written as 

h= y, Jlsisn i+1 + jzY^-s^+ H Y. s ?’ ( S4 ) 

i,~y=x,y ,z i i 


where S 7 is the 7 component of Sj. The xy components of the first term can be divided into an exchange term 


JT+J? 


(S+S~+h.c.) and a double spin-flip term 




( 5 , 


+ o+ 
‘-’i+l 


4 —^ i+1 + h.c.). This term has the same form as the order 
). It might be dominant at high magnetization M/M s ~ 1. Therefore, 


4 V'-'i ^ i -\-1 

parameter of the nematic state (£7 ~Sr + 1 
qualitatively a nematic state might be expected near the saturation in the presence of an xyz-exchange anisotropy. 
More details will be published elsewhere. 
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